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I. INTRODUCTION 


This final report gives details of the first year 
of a four year numerical ocean circulation modeling program 
for the Gulf of Mexico. The aim of the program is to pro- 
gressively upgrade, in modest increments, an existing numer- 
ical ocean circulation model of the Gulf so that the final 
model has a horizontal resolution of about 10 km and verti- 
cal resolution approaching 1 to 10 m in the mixed layer, 10 
mat the thermocline and 100 m in the deep water. Through- 
out the four year period, the validity of the upgraded model 
will be continuously tested, and velocity field time series 
delivered periodically based on the most realistic simula- 
tion of Gulf circulation available. 

Experiments in the first year were with the exist- 
ing NORDA/JAYCOR two layer hydrodynamic primitive equation 
ocean circulation model of the Gulf ona 0.2 degree grid. 
They concentrated on correctly specifying the cuastline and 
bottom topography for maximum realism in circulation simula- 
tion, and on how best to include wind forcing. Details of 
selected experiments are presented in this report. 

Experiment 9 represents the best (compared to our 
incomplete knowledge of the real Gulf) simulation available 
at the beginning of the project. It is forced by flow 
through the Yucatan Straits only (no wind forcing), and 
exhibits many of the flow features observed in the Gulf. 
Simulated surface currents sampled every ten days for three 
Loop Current eddy cycles (1140 days) were delivered to MMS 
at the start of the contract period as an ‘early simulation 


run'. Experiment 34 is similar to Experiment 9, but with 


the addition of seasonal wind forcing. The basic circula 

tion patterns show far more variability in this case. Ex 

periment 40 has no wind forcing and its total inflow trans 

port is identical to that in Experiment 9, but the distribu 

tion of transport between the model's two layers has been 
changed (upper layer transport reduced). It exhibits Loop 
Current eddies that are nearer to the size observed in the 
Gulf (Experiment 9 has rather large eddies). The irer<cased 
lower layer flow helps prevent intrusion onto shelf areas, 
and its sea surface variability is remarkably similar to 
that obtained from satellite altimeter crossovers for the 
Gulf. Experiment 60 is identical to Experiment 40 except 
that the horizontal eddy viscosity has been reduced. Some 
of the flow features seen in Experiment 9 were less obdDviosus 
in Experiment 40, but the latter's lower velocities allowed 
the reduction in eddy viscosity and Experiment 60 exhibits 
these features plus some new circulation patterns. Experi- 
ment 68 ic identical to Experiment 40 with the addition of 
wind forcing from the Navy Corrected Geostrophic Wind data 
set for the Gulf. This wind set has wind stresses every i2 
hours from 1967 to 1982. The addition of winds increases 
the velocities encountered, and attempts to add this wind 
forcing to Experiment 60 were not successful. Simulated 
surface currents sampled every three days for more than 10 
years were delivered to MMS from Experiment 68, as repre- 
senting the best simulation available from the first year 
effort. 


II. EXPERIMENT 9 


Since the ocean model contained many innovative 
features it was discussed in detail in Hurlburt and Thompson 
(1980). In particular Section 2 (pp 1613-1614) gives the 
model equations and Appendix B (pp 1647-1650) describes the 
numerical design of the model. Since that time the capabil- 
ity to handle general basin geometry has been added but this 
does not affect the description in any major way. Wind 
forcing is treated identically to interfacial and bottom 
Stress terms, i.e., wind stress appears directly as an addi- 
tive term in the momentum equation [see p 1614 of Hurlburt 
and Thompson (1980)]. 

In terms of ‘realism’ Experiment 9 was the most 
successful Gulf of Mexico numerical simulation prior to the 
Start of this project. The model was driven from rest to 
Statistical equilibrium solely by a steady inflow through 
the Yucatan Straits which was compensated by outflow through 
the Florida Straits. The model parameters were: 


O upper/lower layer inflow transport = 26/4 Sv, 

© horizontal eddy viscosity, A = 300 m*/sec, 

o Coriolis parameter at the southern boundary f = 5 x 
10°? sec’, 

re) gravitational acceleration, g = 9.8 m/sec’, 

© reduced gravity, g' = .03(H1+H2)/H2 m/sec’, 

o reference layer thicknesses, Hi = 200m and H2 = 
3400m, 

Oo minimum depth of bottom topography = 500m, grid 


spacing, 25 by 25 Km, 


beta, df/dy = 2x 107'' m' sec’, 


wind stress = 0O, 
interfacial stress = 0, 
coefficient of quadratic bottom stress = .002; and 


o oOo 8 98 °@O 


time step = 1 hour. 


Figure 1 compares ‘instantaneous’ upper ocean flow 
patterns just before an eddy is shed by the Loop Current (a) 
from the numerical model and (bd) from observations by Leip- 
per (1970). The ability of the model to simulate observed 
features is clearly demonstrated by this comparison, which 
is remarkable given the simplicity of the model forcing. 
However some discrepancies remain, for example the eddy has 
not penetrated as far into the Gulf and is more intense than 
that shown in the observations. Waves can be seen moving 
around the wall of the Loop Current in both the model and 
the observations, but in the model they are at the limit of 
resolution and erefore unrealistically large. Moreover in 
the Gulf the waves are more pronounced on the eastern wall 
of the Leop and can form strong cold intrusions that may 
contribute to the eddy shedding process (Vukovitch and Maul, 
1984). This is an example of a feature that would benefit 
greatly from 10km model grid resolution. As shed eddies 
propagate westward (Figure 2a) the model spontaneously de- 
velops a counter rotating vortex pair (Figure 2b), a struc- 
ture repeatedly observed in the Western Gulf (Figure 3). 
The roles of the wind and the Loop Current eddies in the 
formation of this structure have been 4 matter of some con- 
troversy (Merrell and Morrison, 1981). Although wind fore- 
ing was not present in this simulation a major role for 


Winds has not yet been ruled out. After spin up the 


experiment sheds an eddy once every 390 days and the eddy 


shedding cycles are very similar. 


IIl. EXPERIMENT 34 


Experiment 34 is similar to Experiment 9, but with 
the addition of wind forcing based on a seasonal climatology 
from ship observations (Elliot, 1979). Linear interpolation 
in time was used between the seasonal fields to produce the 


wind stress at each time step. The model parameters were: 


© upper/lower layer inflow transport = 26/4 Sv, 

© horizontal eddy viscosity, A = 300 m*/sec, 

© Coriolis parameter at the southern boundary f = 4.5 
x 10°? sec’, 

© gravitational acceleration, g = 9.8 m/sec’, 
reduced gravity, g' = .03(H1+H2)/H2 m/sec’, 
reference layer thicknesses, Hi =» 200 mand H2 «= 


3300 m, 


© minimum depth of bottom topography = 500m, 

© grid spacing, 20 by 22 Km (0.2 by 0.2 degrees), 

o beta, df/dy = 2x 107'' m' sec’, 

o wind stress from seasonal climatology based on ship 
observations, 

© interfacial stress «= 0, 


coefficient of quadratic bottom stress = .002; and 
o time step = 1.5 hours. 


Seasonal climatological winds obdviously cannot 
represent most of the wind variability, and Experiment 31 
which was driven solely by these winds attains a steady 
yearly cycle (JAYCOR, 1983). Even s0 the addition of wind 
forcing increases the variabiiity of the Loop Current 


system, including the eddy shedding period and eddy path. 
For examxle Figure 4 compares 360 model days from Experi- 
men:s 28 and 34 (which are identical except that Experiment 
26 has no wind forcing). From these snapshots, taken every 
90 days, there is little difference between the two experi- 
ments. But if Experiment 34 is sampled ever 20 days, as in 
Figure 5, it is apparent that eddies were shed in the space 
of about ont ear. Figure 5 also shows that the circulation 
pattern in the western Gulf can change very rapidly at 
times. Figurzw 6 shows the mean interface deviation and its 
variability for experiments with wind forcing only, port 
forcing only and with wind plus port forcing. This demon- 
Strates that even in the mean the interaction of wind and 
port forcing is not linear, i.e., the mean of the dual fore- 
ing experiment is not the sum of the other two means. The 
variability is increased in the dual forcing case, particu- 
larly in the central and western Gulf. 

Figures 4, 5 and 6 also clearly demonstrate that 
ocean circulation climatologies are inappropriate for use in 
Oil spill risk analysis in the Gulf of Mexico. Figures 4 
and 5 show the highly dynamic mature of fiow in the Gulf, 
and that there is no obvious averaging period (since the 
dominant Loop Current eddy shedding cycle can take anywhere 
from 6 to 18 months). Figure 6 shows that variability is as 
strong as the mean signal, and that the mean is quite dif- 


ferent from circulation at any given point in time. 


IV. EXPERIMENT 40 


Experiment 9 exhibits many of the circulation fea- 
tures found in the Gulf, but discrepancies remain. Some 
problems, such as the correct simulation of the waves on the 
wall of the Loop Current, can only be solved by upgrades to 
the model proposed for years 2, 3 and 4of this project. 
Others, such as the correct simulation of circulation on 
shelf areas, cannot be completely solved without a break- 
through in modei design which is outside the original scope 
of the project. However the major aim of the first year 
effort was to investigate just how realistic the simulation 
could be made without major changes to the model. 

Loop Current intrusion onto the continental shelf 
was identified as a major problem area. This is caused by 
the fact that the model's bottom topography is confined to 
the lower layer, so the minimum topography depth is taken to 
be 500 m and there is a flat shelf in the model between the 
500 m isobath and the coast (Figure 7). Note that intrusion 
of strong Loop Current related flows onto the continental 
slope/shelf «.es occur in the Gulf, so the problem is how to 
control such intrusions in the model and how to have confir 
dence in the results given the apparent deficiencies in 
model formulation. 

Coastai areas are so important for oil spill risk 
analysis that the use of a layered circulation model might 
be carefully examined. But the Loop Current and its associ- 
ated eddies dominate the overall Gulf circulation, and have 
amajor impact on shelf circulation in both the eastern and 
western Gulf. Given the state of the art in supercomputers 


it is simply not practical to produce a circulation model of 
the Gulf that has 10 Km (or even 25 Km) horizontal grid 
resolution Gulf-wide to simulate the Loop Current systen, 
and high vertical resolution to give improved simulations 
over the shelf. Therefore the choice was between a layered 
model such as the one used here (possibly coupled to local 
shelf models), or a level model not significantly better 
Gulfewide than the existing geostrophic climatology. In 
year four of the project, the layered model will be coupled 
to a one dimensional mixed layer model with high vertical 
resolution. This will improve simulation accuracy over 
shelf areas, but will not solve the intersection problem. 
Further details on the question of model design can be found 
in the original proposal (JAYCOR, 1983). 

Several ad hoc methods were tried to control the 
flow over the shelf areas. For example the Yucatan coaste 
line was extended te cover most of the Campeche Bank in an 
attempt to prevent early westward bending of the Loop Curs* 
rent, and interfacial friction was applied over shelf areas 
Only to try to control intrusion. But none of these at* 
tempts were very successful, and the best solution came from 
Simulations addressing the fact that the Loop Current eddies 
in Experiment 9 are large and have high maximum currents 
(they are at the very edge of the acceptable range of eddy 
sizes). 

Eddy radius is dependent on the upper layer velocir- 
ty at the core of the Loop Current (Hurlburt and Thompson, 
1980). It can be controlled by three parameters, (a) upper 
layer rest thickness, (b) the density contrast between lay~ 
ers (i.e., g@'), and (¢c) inflow transport and its distribu* 


tion between layers. Upper layer thickness and g' were 


carefully chosen in the original experiments to give the 
best representation possible of the Gulf in a two layer 
hydroauynamic model. Data on the actual inflow transport 
through the Yucatan Straits is not plentiful, but the figure 
of 30 Sv for the total average transport is consistent with 
what data is available. Data on the distribution of that 
transport in the vertical is almost non-existent, indeed 
even the direction of deep flow is not entirely certain. 
Therefore the original aistribution of 26 Sv in the upper 
layer and 4 Sv in the lower layer was somewhat arbitrary, 
and the upper layer transport can be lowered to produce 
Smaller eddies. Exactly what range of eddy sizes is realis= 
tic is hard to quantify, but there is one source of Gulf 
wide data that can be used as a guide. Maps of sea surface 
variability for the Gulf nave been produced from all hydro- 
graphic, STD and XBT data (Maul and Herman, 1984), and from 
Satellite altimeter cross-overs (Marsh, et al., 1984). The 
20 Sv upper and i0 Sv lower layer distribution of inflow 
transport in Experiment 40 gives rise to a variability map 
very similar to that obtained from the satellite (Figure 8), 
these maps agree more closely with each other than with the 
map from hydrographic data (Figure 9). Based on the agree=- 
ment of variability maps, the mean sea surface from Experi- 
ment 40 may well be the best mean available for the Gulf 


(Figure 10). The model parameters were: 


© upper/lower layer inflow transport = 20/10 Sv, 
© horizontal eddy viscosity, A = 300 m*/sec, 
© Coriolis parameter at the southern boundary f = 4.5 
35 al | 
x 10 sec , 


© gravitational acceleration, g = 9.8 m/sec’, 


10 


© reduced gravity, g' =°.03(H1+H2)/H2 m/sec’, 
o reference layer thicknesses, Hi = 200 m and H2 = 
3300 o, 


minimum depth of bottom topography = 500 o, 
grid spacing, 20 by 22 Km (0.2 by 0.2 degrees), 
beta, df/dy = 2 x 10°!! n”! sec”!, 


wind stress = Q, 
interfacial stress = 0, 


coefficient of quadratic bottom stress = .003; and 


o 9° 8 8 8 8 °@O 


time step = 1.5 hours, 


As a side effect of increased flow in the lower 
layer intrusion onto the Florida Shelf has been reduced, as 
can be seen from a comparison of free surface snapshots from 
Experiments 28 (with 26/4 transport distribution) and 40 
(Figure 11). This is because lower layer flow tends to fol- 
low the bottom topography contours (or more exactly f/H con- 
tours), and the increased deep transport allowed the upper 


layer flow to feel the continental slope more strongly. 


a 


V. EXPERIMENT 60 


Horizontal eddy viscosity is used in the model to 
parameterize subegrid scale processes. As a general rule 
eddy viscosity should be chosen as 1oOw as possible in high 
resolution models, although if it is too low circulation 
features can be produced that are not adequately resolved by 
the model grid. These features will not necessarily be 
Simulated accurately. For example, the waves moving around 
the Loop Current eddy in Experiment 9 fall into this cate 
gory. There is no substitute for high horizontal resolue 
tion. 

Experiment 40 does not exhibit the smaller sca]? 
features, such as the meanders on the wall of the Loop Cur- 
rent and the counter-rotating vortex pair as dramatically as 
earlier experiments. However since this experiment has 
lower maximum speeds the horizontal eddy viscosity can be 
lowered. Experiment 60 therefore is identical to Experiment 
40 except that the eddy viscosity is 100 rather than 300 


m*/sec. The full model parameters were: 


Oo upper/lower layer inflow transport = 20/10 Sv, 

© horizontal eddy viscosity, A = 100 m*/sec, 

© Coriolis parameter at the southern boundary f = 4.5 
x 10°? sec”’, 

© gravitational acceleration, g = 9.8 m/sec’, 

© reduced gravity, g' = .03(H1+H2)/H2 m/sec’, 

o reference layer thicknesses, Hi = 200 m and H2 = 
3300 m, 

Oo minimum depth of bottom topography = 500 m, 


12 


grid spacing, 20 by 22 Km (0.2 by 0.2 degrees), 


beta, df/dy = 2 x 107! a”! sec” 


wind stress = 0, 


interfacial stress = QO, 


coefficient of quadratic bottom stress = .003; and 


o oOo 8 8 8 98 


time step = 1.5 hours. 


Figures 12 to 19 show upper layer velocities from 
Experiment CO covering 300 days at irregular intervals. 
Figures 12, 13 and 14 are 40 days apart and show a Loop 
Current eddy moving into the western Gulf. In Figure i2 
there is a pair of eddies in the northeast Gulf. In Figure 
13 there are six or more small eddies in the northeastern 
Gulf that have been spun off the main Loop Current eddy, 
this may be the feature observed as a meandering current in 
Brookes and Legeckis (1984). Forty days later the currents 
in this area have again changed and most of the flow is to 
the east (Figure 14). Figure 15 shows the situation one 
hundred days later. The flow patterns in the entire western 
half of the Gulf are extremely complex as the Loop Current 
eddy dissipates on the coast of Mexico at about 24N. Note 
also the current along the continental slope off Florida. 
Ninety days later an eddy is about to break off from the 
Loop Current, Figures 16 to 19 are ten days apart. Cyclonic 
rings are moving around the wall of the Loop, and one ap* 
pears to pass right through the space between the Loop and 
the shed eddy (Figure 19). During this time the eddy is 
intruding significantly onto the Campeche Bank. 

On 19 November 1980, the NOAA Buoy Office (NDBO) 
deployed three experimental TzD drifting buoys, at approxi- 
mately 24.5N and 92W, in an eddy that had just been shed by 


13 


the Loop Current (Kirwan et al., 1984). All three drifters 
stayed in the eddy for at least five months as it propagated 
westward to the Mexican coast, see Figure 20. The buoys 
were undrogued but had 200 m thermistor cables, which clear# 
ly coupled the drifters to the deeper circulation. Figure 
21 shows simulated drifter tracks for 160 days from Experi- 
ment 60. They start earlier in the eddy cycle than the NBDO 
buoy tracks. The simulated drifter woves in response to the 
upper layer velocity from the ocean model, which represents 
the mean velocity above the thermocline. Along the drifter 
tracks the upper layer thickness is between 250 and 350 om. 
The simulated drifter position is calculated every 45 min# 
utes based on the velocities at the four previous positions 
using a AdamseBashforth prediction method. Velocities are 
linearly interpolated in space and time between the archived 
model velocity fields that are available once every ten 
model days. Interpolation in time is not an ideal way to 
Calculate velocities since the eddy moves west about 30 Km 
in ten days. The slight elongation of the eddy's east-west 
axis this causes is acceptable, but it also filters out any 
short time scale velocity fluctuations. One advantage of 
Simulated drifters is that the ocean model data also gives a 
view of the entire Gulf, Figures 22 and 23 show snapshots of 
the model free surface deviation every 30 days from the 
simulated deployment date (model day 1680) until model day 
2010. 

The simulated drifters are in good general agree-~ 
ment with the actual drifter data. Both follow approximate- 
ly the same path into the southwest Gulf. Experience with 
the model suggest that this is the preferred eddy trajectory 


although they can also track due west to arrive at the coast 
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further north. The observed average rotation period is 
between 14 and 17 days, with a westward translation speed of 
5 to 10 cm/sec and velocity component speeds of on the order 
of 50 em/sec (Kirwan et al., 1984). The simulated eddy has 
a rotation period of 15 to 16 days, a westward translation 
speed of 3 to 6 em/sec and velocity component speeds of on 
the order of 50 cm/sec. The lower model eddy translation 
speeds may be due to the absence of wind forcing, or to the 
model's idealized vertical structure. The addition of a 
third layer and thermodynamics to the model in year three of 
this project will improve the vertical density structure, 
and may lead to different translation speeds. 

The simulated drifters exhibit very regular loops 
and appear to remain at approximately the same relative 
position within the eddy for long periods. The actual 
drifters on the other hand are much more variable, with 
paths indicative of changes of the drifter location relative 
to the ring center. These changes are probably primarily 
due to windage effects on the drifters, and might be mini- 
mized by adding a drogue at the end of the thermistor line 
(i.e., at 200 m) on future buoys. In principle windage 
could be accounted for in drifter simulations; model experi-« 
ments with wind forcing automatically account for layer 
averaged Ekman effects, but allowing for winds acting di« 
rectly on the buoy would require knowledge of the the rela- 
tive effectiveness of such forcing. 

Figure 24 follows the path of simulated drifter 
number 3 for 300 days. Remarkably it is still tracking the 
eddy after all this time, even though the eddy has almost 
totally dissipated by model day 1980 (see Figure 23). 
Drifter 1599 tracked the eddy remnants until mid-June 1981, 


15 


but drifters 1598 and 1600 left the eddy in mid-April and 
early=May respectively. The model eddy probably dissipated 
too slowly, because its interaction with the continental 
shelf cannot be modeled accurately. But equally it is not 
necessarily the case that the actual eddy had entirely dis- 
appeared by mid*June 1981, it is possible that wind effects 
caused drifter 1599 to leave the eddy at that time. In any 
case the model accurately simulates the northward motion of 
the eddy once it reaches the coast of Mexico. However the 
paths of 1598 and 1600 once they leave the eddy suggest that 
the remnant of a previous Loop Current eddy that persists 
off the Texas coast in the model simulation throughout this 
time period was not present in the summer of 1981. Similar 
features occur in almost every model simulation, even in 
Simulations with wind forcing only, where it is a wind ine 
duced gyre rather than a Loop Current eddy. Their presence 
is explainable by the northward migration of anti*cyclonic 
rings along the coast until they reach Texas, where the 
continental slope turns east and they can go no further. 
But in the Gulf the rings probably dissipate more quickly 
against the continental shelf than they do in the model, 
Since the latter cannot include topography shallower than 
500 m. 

Figure 25 contains time series of velocity for 
drifters 1598, 1599 and 1600. The strong high frequency 
contribution is unusual and is largely due to the 30 hour 
basin tidal resonance and possibly inertial oscillations, 
along with diurnal and semidiurnal tides and a 7-hour free 
gravity mode (Kirwan et al., 1984). Figure 26 shows time 
series of velocity for simulated drifters 3 and 4. These 


show no high frequency components, as is to be expected 


16 


given that the simulation only has access to new model velo-=- 
city fields once every ten days. High frequency components 
might appear in simulated drifters that were calculated "on 
the fly" within the ocean model, since new velocity fields 
would then be available every 90 minutes. The low frequency 
velocity components of the actual buoys agree well with the 
simulated drifters. Both show periodicity associated with 
the etdy circulation and velocities of about 50 cm/sec. The 
Simulated time series are far more regular, as is expected 
from the comparison of drifter tracks. 


17 


Vi. EXPERIMENT 68 


Simulations forced by winds based on a seasonal 
Climatology from ship observations have already been de- 
scribed. Such wind fields are not ideal for driving ocean 
models since they contain very little of the total wind 
variability and mean wind strengths are in general far weak= 
er than instantaneous winds. Recognizing this deficiency 
NORDA funded JAYCOR to produce a wind set for the Gulf 
based on the Navy's twelve hourly surface pressure analysis, 
which is available from 1967 to 1982 (Rhodes et al., 1984). 
The geostrophic winds, corrected geostrophic winds, and wind 
stresses (all on a one degree grid covering the Gulf) every 
12 hours from 1967 to 1982 are on magnetic tape. These will 
be made available through the MMS Gulf of Mexico regional 
office. 

Figure 27 shows the wind stress and wind stress 
curl from this data set for 0000 and 1200 GMT on 14 January 
and 0000 GMT on 15 January 1976. There is large temporal 
variability of the wind field during this period, as general 
easterly flow gives way to strong northerly flow after a 
frontal passage in just a 24 hour period. The wind stress 
curl field also shows the rapid change, from a relatively 
weak field to a very strong field with strong horizontal 
gradients. Figure 28 shows similar plots for 14 and 15 July 
1976. Even in the summer, when flow is generally weaker, 
very significant differences can be seen in a short time 
period. These strong variations and very rapid changes in 
the wind field indicate why the modeling of Gulf circulation 


requires wind data on short temporal scales. 
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Figures 29 to 32 show the seasonal climatologies 
averaged over the period 196741982. The wind stress and 
wind stress curls are much stronger in the winter sez2son 
than the summer season as would be expected. There are 
persistent areas of positive curl over the Yucatan and nega~ 
tive curl in the southwest Gulf that are present for all 
seasons, but were not seen in any previous study of Gulf 
wind stresses. Although not present at all time periods 
(Figures 27 and 28), these are also the dominant features of 
the instantaneous curl fields. 

There have been no published accounts of driving an 
ocean circulation model for long time periods with winds 
sampled as frequently as those available in this wind set. 
Previously a monthly climatology would have been considered 
an exceptionally good data set for such an application. One 
of the initial goals of this first year effort was to deter- 
mine how best to use this data set to drive the ocean model. 
The model simulates the layer averaged circulation above and 
below the thermocline, and so only includes the longer term 
effects of winds on ocean currents. But it is clearly 
preferable for the ocean model to integrate the effect of 
short term variability on these currents, rather than for 
the winds to be averaged before input to the model. Before 
testing began it was expected that inertial oscillations and 
gravity waves generated by the highly variable forcing would 
make using the 12 hourly winds directly impractical. So the 
plan was to test the model with increasingly long wind aver~- 
ages until these problems became manageable. In all cases 
the wind stresses are linearly interpolated in time between 
inputs, so the wind forcing is slightly different at each 
time step. 


19 


Experiment 68 is identical to Experiment 40 except 
for the addition of wind forcing after the port forced cir* 
culation has fully spun up. Wind input is every i2 hours, 
at first 1967 winds were used repeatedly to spin up the 
wind driven flow and then winds from 1967 to 1977 were ap- 
plied in sequence. The expected difficulties with frequent 
wind input did not arise, although attempts to add these 
winds to Experiment 60 (with lower eddy viscosity) were 


unsuccessful. The model parameters were: 


upper/lower layer inflow transport «= 20/10 Sv, 
horizontal eddy viscosity, A = 300 m*/sec, 
© Coriolis parameter at the southern boundary f = 4.5 
m5 » 
x 10 sec, 
gravitational acceleration, g = 9.8 m/sec’, 


reduced gravity, g' = .03(H1+H2)/H2 m/sec’, 


Oo reference layer thicknesses, Hi = 200 mand H2 = 
3300 mo, . 

o minimum depth of bottom topography = 500 o, 

© grid spacing, 20 by 22 Km (0.2 by 0.2 degrees), 

Oo beta, df/dy = 2 x 107!! a”! sec”, 

o wind stress from 12 hourly Navy Corrected Geostror- 
phic Wind set, _ 

o interfacial stress = 0, 


coefficient of quadratic bottom stress = .003; and 
o time step « 1.5 hours. 


Figures 33 to 39 show upper layer currents (i.e., 
vertically averaged currents above the thermocline) every 60 
days for 360 days. Vectors are only drawn at every second 
point (i.e., every 0.4 degrees) to improve readability. 
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Figure 34 shows the furthest northward penetration of the 
Loop Current ever attained by the ocean model, this configu* 
ration is often seen in the Gulf. After the eddy breaks off 
the Loop Current intrudes onto the Florida Shelf and some of 
the flow splits off to the north for a brief time (Figures 
35 and 36). Similar intrusions have been observed in the 
Gulf, but the models inadequate representation of shelf 
topography make it likely that the simulated currents in 
shallow areas (say less than 100 m) are too high. A persis-~- 
tent antiecyclonic gyre in the north west Gulf has been a 
feature of almost all Gulf simulations performed to date. 
The additicn of wind forcing in Experiment 68 has increased 
its average size and its effect on incoming Loop Current 
eddies (Figures 37, 38 and 39). The presence of a gyre in 
this position is explainable by the northward migration of 
antiecyclonic eddies along the coast of Mexico until the 
continental slope bends eastward and they can go no further. 
However in the Gulf the gyre probably dissipates relatively 
rapidly against the shallow shelf area. The ocean model 
cannot include a shallow shelf, but it may be possible for 
it to parameterize the effect of such a region on strong 
currents. Ad hoc patches to the model of this kind can only 
be justified if (a) there is clearly a problem with the 
existing model simulation and (b) there is sufficient obser- 
vational data available to verify that the patch is indeed 
"correcting" the simulation towards more realistic flow 
patterns. Modifying the simulation towards the general 
perception of what the circulation should be, without sup- 
porting data, is dangerous. Our understanding of Gulf cir- 
Culation has changed radically in the last ten years, and 


given the sparsity of the existing observational data base 
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it is likely to continue changing. In the case of the gyre 
in the northwest Gulf condition (a) probably can be satisfir 
ed, but obtaining the data to satisfy condition (b) will be 
difficult. MMS's planied observational program in the west@ 
ern Gulf will be of assistance in this area. Finally it 
should be noted that problems of this kind can sometimes be 
solved by less drastic measures. For example the Loop Cur- 
rent intrusions onto the Florida Shelf were eventually ree 
duced by re-evaluating the inflow transport distribution. 
Similarly the north western gyre was least odvious in Exper- 
iment 60 which has the lowest eddy viscosity. It is possi- 
ble therefore that next years simulations with a 10 km grid 
and correspondingly lower eddy viscosity will resolve this 
question. 

Simulated surface currents from Experiment 68 have 
been delivered to MMS as representing the best simulation 
Gata available to date. They consist of velocity component 
(u and v) fields on a 0.2 degree rectangular grid covering 
the Gulf area, sampled every three days for 3780 model days 
(10.3 years). Velocity plots taken once every 30 days (i.e. 
plots of every tenth set of fields) for the entire ten year 
period were also delivered. Figures 33 to 39 represent a 
very small subset of these plots, note that a current vector 
is only drawn at every second point (i.e., on a-0.4 degree 
grid) to improve readability. 

The velocities represent wean currents above the 
thermocline and are therefore essentially independent of 
local winds. To a good approximation they can be treated as 
geostrophic surface currents. The Oil Spill Trajectory 
Analysis (OSTA) model can therefore include the simulated 


Surface currents in exactly the same way as the previous 
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Climatologica) geostrophic surface currents. Except that 
now new current data jis available more frequently, and the 
OSTA model can include variation in surface turrents (over a 
ten year period) as well as variation in wind driving. 
However the use of ten years of surface currents, rather 
than climatology (which is effectively one year of data), 
means that the OSTA model will require more computer re- 
sources to complete its statistical analysis of risk. 
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FIGURE 1: 


FIGURES 1-39 


(a) Instantaneous view of the interface deviation 
in a two-layer simulation of the Gulf of Mexico 
driven from rest to Statistical equilibrium 
solely by inflow through the Yucatan Straits 
(Experiment 9). The contor interval is 25 o, 
with solid contours representing downward devia- 
tions. (b) Depth of the 22 degree isothermal 
surface, 4-18 August 1966 (Alaminos cruise 66-A- 
11), from Leipper (1970). The contour interval 
is 25 m. 
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FIGURE 2: 


(a) Interface deviation from the Gulf of Mexico 
Simulation at model day 1970 after an eddy has 
separated from the Loop Current and propagated 
westward. (b) Ninety days later the major anti- 
cyclonic eddy at day 1970 has developed into a 
counter-rotating vortex pair in the western Gulf. 
The cyclonic vortex is to the north and the anti- 
cyclonic to the south. 
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FIGURE 3: 


Counter-rotating vortex pair in the western Gulf 
of Mexico as shown by the depth of the 15 degree 
isotherm (in meters), observed in April 1978. 
The cyclonic vortex is to the north and the anti- 
cyclonic to the south (from Merrell and Morrison, 
1981). 
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FIGURE 4: 


Instantaneous view of the interface deviation 


every 90 days, from day 90 of model 
0 of model year 10, for Experiment 
Experiment 34 (right). Experiment 
cal to 28 except for the addition 
ing. The contour interval is 25 m. 


year 9 to day 
28 (left) and 
34 is identi- 
of wind force 
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FIGURE 5: Instantaneous view of the interface deviation 
every 20 days, from day 260 of model year 9 to 
day 9 of model year 10, for Experiment 34. 
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FIGURE 6: 


Interface deviation mean and variability, for the 
Gulf of Mexico from ocean model. (a) Experiment 
31, wind forcing only; (b) Experiment 28, port 
forcing only; (c) Experiment 34 wind plus port 
forcing. The contour interval is 12.5 m. 
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FIGURE 7: 


Bottom topography and coastline geometry for Gulf 
of Mexico model on 0.2 degree grid. The contour 
interval is 250 m and the shallowest depth is 500 
m. The section of the Caribbean shown is treated 
as land by the model, the position of the inflow 
port is marked by the termination of contour 
lines in the Yucatan Strait. 
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FIGURE 8: 


Sea surface height variability for the Gulf of 
Mexico. (a) Based on about 16,000 GEOS-3 and 
SEASAT cross overs, spanning nearly four years 
(from Marsh, Cheney and McCarthy, 1984). (b) 
Based on an ocean model simulation with port 
forcing only (Experiment 40), measured over three 
eddy cycles at statistical equilibrium with the 
free surface sampled every ten days for a total 
of over 300,000 "observations". 
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SEA SURFACE VARIABILITY FROM 
GEOS-3 AND SEASAT CROSS OVERS (CM) 
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FIGURE 9: Sea surface height (a) variability and (b) mean, 
for the Gulf of Mexico. Based on all available 
hydrographic, STD and XBT data at over 16,000 


Stations, with substantial filtering (from Maul 
and Herman, 1984). 
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FIGURE 10: 


Mean sea surface height for the Gulf of Mexico. 
Based on an ocean model simulation with port 
forcing only (Experiment 40), measured over 
three eddy cycles at statistical equilibrium. 
The contour interval is 5 cm. 
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FIGURE 11: 


Instantaneous view of the free surface deviation 
(a) from Experiment 28 (with 26 upper and 4 Sv 
lower layer inflow transport), and (b) Experi- 
ment 40 (with 20 upper and 10 Sv lower layer 


inflow transport). In similar phase of eddy 
cycle, Experiment 40 shows less intrusion of the 
Loop Current onto the Florida Shelf. The 


contour interval is 10 cm. 
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FIGURE 12: Instantaneous view of upper layer averaged velo- 
cities from Experiment 60 on model day 2130, 
velocities above 50 cm/sec are not shown. 
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FIGURE 13: Instantaneous view of upper layer averaged velo- 
cities from Experiment 60 on model day 2170, 
velocities above 50 cm/sec are not shown. 
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FIGURE 14: Instantaneous view of upper layer averaged velo- 
cities from Experiment 60 on model day 2210, 
velocities above 50 cm/sec are not shown. 
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GE@STR. CURRENTS G. OF MEXICO 0, 60 
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FIGURE 15: Instantaneous view of upper layer averaged velo- 
cities from Experiment 60 on model day 2320, 
velocities above 50 cm/sec are not shown. 
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FIGURE 16: 


Instantaneous view of upper layer averaged velo- 
cities from Experiment 60 on model day 2410, 
velocities above 50 cm/sec are not shown. 
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FIGURE 17: Inastantaneous view of upper layer averaged velo-~ 
cities from Experiment 60 on model day 2420, 
velocities above 50 cm/sec are not shown. 
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GE@STR. CURRENTS G. OF MEXICO 0. 60 
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FIGURE 18: Instantaneous view of upper layer averaged velo- 
cities from Experiment 60 on model day 2430, 
velocities above 50 cm/sec are not shown. 
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FIGURE 19: Instantaneous view of upper layer averaged velo- 
cities from Experiment 60 on model day 2440, 
velocities above 50 cm/sec are not shown. 
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FIGURE 20: 


Paths of NDBO drifters (a) 1598, (b) 1599, and 
(c) 1600 from November 20, 1980, through May 11, 
1981. The numbers 0 through 6 give the posi- 
tions on November 20, December 20, January 20, 
February 20, March 20, April 20, May 11, respec- 
tively (from Kirwan et al., 1984). 
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FIGURE 21: 


Paths of simulated drifters (a) 3, ané (bd) 4 
from model day 1680 to model day 1840 of Gulf 
model Experiment 60 (which is forced solely by 
inflow through the Yucatan Straits). The tracks 
Start earlier in the eddy cycle than those in 
Figure 20. The simulated drifter moves in re- 
sponse to the upper layer velocity from the 
ocean model, which represents the mean velocity 
above the thermocline. Along the drifter tracks 
the upper layer thickness is between 250 and 350 
m. The track is drawn as a solid line for 20 
days, then dashed for 20 days, and so on. There 
is a dot every 5 days. 
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FIGURE 22: 


Instantaneous view of the free surface deviation 
every 30 days, from model day 1680 to model day 
1820, for model Experiment number 60. The con- 
tour interval is 10 ecm, and solid contours re-~ 
present upward deviations with respect to the 
sea surface height at rest. 
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FIGURE 23: Instantaneous view of the free surface deviation 
every 30 days, from model day 1860 to model day 
2010 (Experiment 60). 
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FIGURE 24: Path of simulated drifter number 3 from model 
day 1680 to model day 1980 (Experiment 60). 
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FIGURE 25: Time series of velocity for drifters (a) 1598, 
(b>) 1599, and (c) 1600. From Kirwan et al., 
(1984). 
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FIGURE 26: Time series of velocity for simulated drifters 
(a) 3 and (bd) 4 from miel day 1680 to model day 
1880 (Experiment 60). 
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FIGURE 27: Instantaneous wind stress and wind stress curl 
from the Navy Corrected Geostrophic Wind data 
set, for 0000 and 1200 GMT on 14 January and 


0000 GMT on 15 January 1976. 
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FIGURE 28: Instantaneous wind stress and wind stress curl 
from the Navy Corrected Geostrophic Wind data 
set, for 0000 and 1200 GMT on 14 July and 0000 


GMT on 15 July 1976. 


78 


WIND STRESS WIND STRESS CURL 


196/1976 AT © GNT ait, 1906/1976 AT © GuT X = 20£-07 
T ss a a ee ee ee ee ee ee 7 7 T T T 3iN 7 a poe | 
‘ ’ PP, . . * _—«_ : 
. ? ‘ “ 
* > * * = — 
4 s . > . > . - - - — =~ » 


4 , A . . . . . . . -_-_—-«< 
~ RA NR BR Be . * « . . . 
ee SL NS @ © «6 © «@ - 


14N * A - 4 ry 4 ry r - s 4 - 
- 985 (DEG) 
MAXIMUM WIND STRESS « 25) DYNES /Cw* 
—amm wv 
WIND STRESS WIND STRESS CURL 
1906/1976 AT i2 GuT at. 1906/1976 AT 12 GuT DC = 20£-07 
316 7 , , - La - 7 AJ A 7 7 7 7 . 7 7 AJ 7 3iN 7 7 7 7 7 , - . 2 7 7 7 7? 7 as 7 7 7 
| “ 
, 8 @© & °°“ ee . » . «os - - 
4 4 a ~ * " . . . ~ —* 
MMA SK 8 6 wy se @ _— 
-_ —<«)-” 


im LoL Pare Paar 
- 985 (DEC) 
MAMMU. WIND STRESS « 3.46 DYNES/cw MAX = 2775-06 
oemme 87) FEC ee wonme 87) (| )-O8r- we 
WIND STRESS WIND STRESS CURL 
1897/1976 AT © GuT aii. 197/1976 AT O GuT DC = 202-07 
ce ee ee ae ce ee oe ee Gee ee 3iN 
VN tte en, . wn ~~ 7 
VA Y «6 [ss wo, =~» 4 
NN ewe cen ame — . 
WQS ARS 
Morsnscxwn ni 
> SS 
SSS SS eo ee 4 (Es) 
SS \ « | 
SSS \ i or 
Lar y Ren es 
be ye NS Ki > Veit: 4 
te ' pf AS 
oo fr imuioesieoteentheh bp 2 4, . -_ 
- ORE (DFG) - 805 OME ~808 


(000. we 


FIGURE 29: Climatological (1967-1982) wind stress and wind 
stress curl from'the ‘Navy Corrected Geostrophic 
Wind data set, for winter (December, January, 
February). 
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FIGURE 30: Climatolcgizal (1967-1982) wind stress and wind 
stress curl from the Navy Corrected Geostrophic 
Wind data set, for spring (March, April, May). 
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FIGURE 31: Climatological (1967-1982) wind stress and wind 
Stress curl from'the Navy Corrected Geostrophic 
Wind data set, for summer (June, July, August). 
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FIGURE 32: Climatological (1967-1982) wind stress and wind 
stress curl from the Navy Corrected Geostrophic 
Wind data set, for fall (September, October, 
November). 
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FIGURE 33: Instantaneous view of upper layer averaged velo- 
cities from Experiment 68 on model day 3858. 
Vectors are only plotted at every second model 
grid point, i.e., every 0.4 degrees. 
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FIGURE 34: Instantaneous view of upper layer averaged veloe- 
cities from Experiment 68 on model day 3918. 
Vectors are only plotted at every second model 


grid point, i.e., every 0.4 gegrees. 
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FIGURE 35: 


Instantaneous view of upper layer averaged velo- 
cities from Experiment 68 on model day 3978. 
Vectors are only plotted at every second model 


grid point, i.e., every 0.4% degrees. 
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FIGURE 36: Instantaneous view of upper layer averaged velo=- 
cities from Experiment 68 on model day 4038. 
Vectors are only plotted at every second model 
grid point, i.e., every 0.4 degrees. 
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FIGURE 37: Instantaneous view of upper layer averaged velo- 
cities from Experiment 68 on model day 4098. 
Vectors are only plotted at every second model 
grid point, i.e., every 0.4 degrees. 
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FIGURE 38: 


Instantaneous view of upper layer averaged velo- 
cities from Experiment 68 on model day 4158. 
Vectors are only plotted at every second model 
grid point, i.e., every 0.4 degrees. 
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FIGURE 39: 


Instantaneous view of upper layer averaged velo 
cities from Experiment 68 on model day 4218. 
Vectors are only plotted at every second model 
grid point, i.e., every 0.4 degrees. 
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APPENDIX A 


MODEL PARAMETERS 


REFERENCE PARAMETERS (EXPERIMENT 40): 


upper layer inflow transport = 20 x 10* m® sec 
(20 Sverdrup), 

lower layer inflow transport = i0 x 10% @ 
(10 Sverdrup), 

wind stress = 0, 

horizontal eddy viscosity, A = 300 m*/sec, 
grid spacing = 20 by 22 km (0.2 by 0.2 degrees), 
upper layer reference thickness, H1 = 200 a, 

lower layer reference thickness, H2 = 3300 a, 
minimum depth of dDsttom tgpography = 500 a, 

beta, df/dy = 2x i0~'’ m' sec’, 

Coriolis pargmeter at the southern boundary, 

f= 4.5 x 10°* sec’', 

gravitational acceleration, g = 9.8 m/sec’, 
reduced gravity, g' = .03 (H1 + H2)/H2 m/sec’, 
interfacial stress = 0, 

coefficient of quadratic bottom stress = .003, and 
time step = 1.5 hours. 


EXPERIMENT 9: 


upper layer inflow transport = 26 x 10° m? sec 
(26 Sverdrup), 

lower layer inflow transport = 4 x 10° m’® sec’', 
(4 Sverdrup), 

grid spacing, 25 by 25 km, 

lower layer reference thickness, H2 = 3400 a, 
Coriolis pgrameter at the southern boundary, 

f= 5 x 107% sec’', 

coefficient of quadratic bottom stress = .002; and 
time step = 1 hour. 

All other parameters as in the reference experi- 
ment. 


EXPERIMENT 28: 


upper layer inflow transport = 26 x 10° m’ sec*! 


(26 Sverdrup), 

lower layer inflow transport = 4 x 10° m’® sec 
(4 Sverdrup), 

wind stress = G, 

coefficient of quadratic bottom stress = .002. 

All other parameters as in the reference experi- 
ment. 


2, 
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EXPERIMcCNT 31: 


@® upper layer inflow transport = 0, 

lower layer inflow transport = 0, 

@e wind stress from seasonal climatology based on ship 
observations, 

e coefficient of quadratic bottom stress = .002. 

@ All other parameters as in the reference experi- 
ment. 


EXPERIMENT 34: 


® upper layer inflow transport = 26 x 10° m® sec 
(26 Sverdrup), 

@ lower layer inflow transport = 4 x 10° m’® sec 
(4 Sverdrup), 

e wind stress from seasonal climatology based on sliip 
observations, 

e coefficient of quadratic bottom stress = .002. 

@ All other parameters as in the reference experi- 
ment. 


EXPERIMENT 40: 

@ All parameters as in the reference experiment. 
EXPERIMENT 60: 

@ horizontal eddy viscosity, A = 100 m*/sec, 

e —" parameters as in the reference experi- 


EXPERIMENT 68: 


@ wind stress from 12 hourly Navy Corrected Geostro- 
phic Wind set, 

@ All other parameters as in the reference experi- 
ment. 
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